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ABSTRACT 

I present a study of the high-energy 7-ray properties of the fiat spectrum radio quasar, 
PKS 1510-089, based on 3.75 years of observations with the Large Area Telescope 
(LAT) detector on-board the Fermi 7-ray Space Telescope. Throughout the observing 
period, the 0.1 GeV < E 1 < 300 GeV 7-ray flux was highly variable, undergoing several 
flaring events where the daily flux exceeded 10 -5 photons cm _2 s _1 on 3 separate 
occasions. The increased photon statistics of these large flares allowed the observations 
to be re-analysed in 6 and 3 hour intervals, revealing flux doubling timescales as small 
as 1.3 ±0.12 hours during the flare rise time, and flux halving timescales of 1.21 ±0.15 
hours during the flare decay. These are the smallest variability timescales measured 
to date at MeV— GeV energies for the Flat Spectrum Quasar class of Active Galactic 
Nuclei. 

The > 10 -5 photons cm _2 s _1 flare events were also studied in more detail in an 
attempt to uncover evidence for the location of PKS 1510-089's 7-ray emission region. 
In particular, two approaches were used: (i) searching for an energy dependence to the 
cooling timescales, and (ii) searching for evidence of a spectral cut-off. The combined 
results of these two approaches, along with the confirmation of > 20 GeV photon 
flux from PKS 1510-089, suggest the presence of multiple 7-ray emission regions being 
located in both the broad line region and molecular torus region of PKS 1510-089. 

An analysis of the highest photon events within the 3.75 year data set finds PKS 
1510-089 to be a source of > 20 GeV 7-rays at the 13.5<r confidence level; a observa- 
tional property which is difficult to explain in the traditional view that 7-ray emission 
from Active Galactic Nuclei originates from the base of the relativistic jet. This gives 
further weight to the argument that there are multiple, simultaneously active 7-ray 
emission regions located along the relativistic jet of Active Galactic Nuclei. 

Key words: galaxies: active - galaxies: individual (PKS 1510-089) - galaxies: jets - 
gamma rays: galaxies. 



1 INTRODUCTION 

The successful launch of the Fermi 7-ray Space Telescope 
affords us an ideal opportunity to investigate the inner work- 
ings of Active Galactic Nuclei (AGN). The ability of the 
Fermi-LAT detector to scan the entire 7-ray sky every three 
hours allows us to study the 7-ray emission from AGN, un- 
biased by activity state or AGN sub-class. This ability has 
revealed the blazar subclass of A GN to be the most numer- 
ous class of known 7-ray sources l|Nolan et all (|2012h h with 
approximately equal number of flat spectrum ra dio quasars 
(FSRQ) and BL Lac objects (|Abdo et al.l (|2011af 0. Further- 
more, not only do blazars dominated the extra-galactic 7- 
ray sky, but during brief periods of intense flare activity, 
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they ca n outshine galactic 7-ray sources as well (|Abdo et al.l 
<2011bltt . 

Apart f rom the radio lobes of the nearby radio galaxy, 
Centarus A (|Abdo et al.l (|2010d )h the 7-ray emission region 
of AGN remains unresolved. As such, the origin of the 7-ray 
emission within blazars remains an open question. Answer- 
ing this question is an active area of AGN research, with 
two main schools of thought: on the one hand, some be- 
lieve the 7-ray emission originates from within the broad 
line region (BLR) of the AGN, while on the other hand, 
some believe the 7-ray originates further out from the cen- 
tral super- massive black hole (SMBH), within the Molecular 
Torus (MT) region of the AGN. 

Traditionally the 7-ray emission has been assumed to 
be close to the base of the relativistic jet, within the BLR. 
The reasoning for this assumption is two-fold: broad-band 
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spectral energy distribution modelling and rapid 7-ray flux 
variability on small timescales. 

Multi-wavelength (MWL) observations of FSRQs have 
found the broad-band spectral energy distribution to be ad- 
equately described by a leptonic model, with the emission 
region located within 300 — 1000 S chwarzschild rad i i from 
the central SMBH (for exam ple, see iGhisellini et "all (|2010h 
and iNalewaiko et al.l (|2012l )'). Likewise, the rapid flux vari- 
ability implies a small emission region size. This is often 
interpreted as evidence of the emission region being located 
close to the base of the jet. This interpretation is based on 
the assumption that the full width of the relativistic jet is re- 
sponsible for the observed 7-ray emission and that the size of 
the emission region, r, is simply related to the opening angle 
of the relativis tic jet, ip, and the d is tance from the SMBH, 
R, via r ~ ipR jDermer et al l {2009), Ghisellini fc Tavecchiol 
l|2009l )1. 

More recently, the wealth of information afford to us 
by the Fermi-LAT detector has found evidence of spectral 
breaks at GeV energ ies in the 7-ray spectrum of some AGN 
l|Abdo et al.l (|2010tJ )). These spectral breaks have been in- 
terpreted in the context of 7-ray absorption through photon- 
photon pair-production w ith the He Lyman r ecomb ination 
continuum of the BLR (|Poutanen fc Sternl (|2010h ). thus 
pointing to a BLR origin of the observed 7-ray flux. It is 
worth noting though, that this interpretation has recently 
been brought into question (|Harris et al.l |2012i 'l'). 

7-ray emission from AGN has also been suggested to 
originate from the MT region of the jet, on the parsec-scale 
distance from the central SMBH. This suggestion is primar- 
ily based on the results of recent MWL observations which 
have found 7-ray flaring events to be accompanied by flare 
events at optical or radio wavelengths, with some of these ra- 
dio flares being resolved o n a parsec-scale distance f r om th e 
SMBH (for example, see Lahtccnmaki & Valtaoka ( 20031). 
iMarscher et"aH (|2010l ). lAgudo et af] (|201ll ), lOrienti et all 

(2orJ). 

Further evidence for an MT origin of the 7-ray emis- 
sion comes from the detection of Very High Energy (VHE) 
7-ray emission from FSRQs. To date, 3 FSRQs have been 
detecte d at energies > 100 GeV, of which PKS 1510-089 
is one JWagnerl ||201Ch . lAleksic et all (|2011al L lATeksic et all 
|2011bl ')'). The photon-rich environment of the BLR of FS- 
RQs is believed to severely attenuate any 7-ray emis- 
sion through photon-photon pair production, resulting in 
a spectral cut-off above > 20 GeV (for example, see 
iDonea fc Protherod (|2003h , lLui fc Bail (|2006h >- As such, the 
detection of VHE emission is difficult to explain with a pure 
BLR-origin for the observed 7-ray. 

One of the primary arguments against an MT origin 
for the 7-ray emission is that the further an emission re- 
gion is from the central SMBH, the bigger it is through 
the combined effects of adiabatic expansion and the open- 
ing angle of the jet. However, this argument is only valid if 
the entire width of the relativistic jet is responsible for the 
observed 7-ray emission. What is more, this argument as- 
sumes that process of adiabatic expansion is dominant over 
any re-collimation process that occurs along the length of 
the jet. The latter assumption is not valid i f the jet un- 
dergoes re-confinement (e.g., see lSokolov et al. (2009)). Fur- 
thermore, detailed computer simulations have found that jet 
instabilities can result in large over-densities in the matter 



distribution at large distance s from t he central S MBH (for 
example, see Nishikawa et al.l (120031') . Perucho et al.l (|2006h . 
iBromberg fc Levinsonl (|2009l ). iKohler et all (|2012l )) 

Interestingly, using the above arguments, evidence has 
been found for the the 7-ray emission region of PKS 1510- 
089 to be located in both the BLR and the MT. From the 
first 11 months of Fermi-LAT operation, Abdo et al. con- 
cluded that^_th«7-rw_emission originated from within the 
BLR (|Abdo et al.l l|2010bh ); this conclusion primarily be- 
ing driven by the presence of a cut-off in the 7-ray spec- 
trum. However, both Marscher et al. and Orienti et al. 
have concluded from their respective MWL campaigns of 
PKS 15 10-089 that the 7-ray o riginates from outsid e BLR 
region (|Marscher et al.l (|2010l 1. lOrienti et all (|2012h 1. with 
their conclusions primarily based on the flaring events being 
observed simultaneously at 7-ray, optical and radio wave- 
lengths. It is important to note however, that these three 
MWL studies focused on three separate flaring events. 

This paper investigates the high-energy 7-ray flux and 
spectral properties of PKS1510-089 during the first 3.75 
years of Fermi-LAT observations. In particular, the in- 
creased photon statistics associated with several large 7- 
ray flare events where the daily flux occasionally exceeds 
10~ 5 photons cm _2 s _1 , allow us to probe flux and spectral 
variability from PKS 1510-089 with unprecedented temporal 
resolution in the MeV-GeV energy range. In §2, I describe 
the Fermi-LAT observations and data analysis routines used 
in this study. The results on flux variability are shown in §3 
with the in-depth flare analysis reported in §4. A brief dis- 
cussion on a multi-zone model is given in §5 with the con- 
clusions given in §7. §6 touches on the VHE 7-ray properties 
of PKS 1510-089 in the context of understanding the origin 
of the observed 7-ray emission. Throughout this paper, a A 
cold dark matter (ACDM) cosmology was adopted, with a 
Hubble constant of Ho = 71 km s _1 Mpc" 1 , Q, m — 0.27 and 
JIa = 0.73 as derived from Wi l kinson Microwave Anisotropy 
Probe results (|Komatsu et alj Hool)). 



2 FERMI-LAT OBSERVATIONS AND DATA 
REDUCTION 

The LAT detector aboard Fermi, described in detail by 
lAtwood et ail (|2009f ). is a pair-conversion telescope, sensi- 
tive to a photon energy range from below 20 MeV to above 
300 GeV. With a large field of view, ~ 2.4 sr, improved an- 
gular resolution, ~ 0.8° at 1 Ge\{]], and large effective area, 
~ 8000 cm 2 on axis at 10 GeV, Fermi-LAT provides an or- 
der of magnitude improvement in performance compared to 
its EGRET predecessor. 

Since 2008 August 4, the vast majority of data taken by 
Fermi has been performed in all-sky-survey mode, whereby 
the Fermi-LAT detector points away from the Earth and 
rocks north and south of its orbital plane. This rocking mo- 
tion, coupled with Fermi-LAT's large effective area, allows 



1 Below 10 GeV photon energy, the 68% containment an- 
gle of the photon direction is approximately given by 9 ~ 
0.8°(E 7 /GeV)~ ' 8 , with the 95% containment angle being less 
than 1.6 times the angle for 68% containment. 
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Figure 1. Lightcurve of the 0.1 GeV < E 1 < 300 GeV 7-ray flux from PKS 1510-089, during the period 2008 August 4 (MJD 54682) to 2012 May 4 
(MJD 56051), binned in daily periods. Only daily periods with a TS > 10 are shown. The light curve of the flare period where the flux exceeded 10 — 5 
photons cm -2 s _1 is shown in the insert. 



Fermi to scan the entire 7-ray sky every t wo orbits, or ap- 
proximately every three hours (|Rit j (|2007l )). This observa- 
tional characteristic of the Fermi satellite affords us, for the 
first time, continuous monitoring of the high-energy 7-ray 
sky, allowing us to study the high-energy properties of AGN 
without suffering the biased of activity state often associated 
with pointed observations. 

The data utilised in this study comprises of all all- 
sky-survey observations taken during the first 3.75 years of 
Fermi-LAT operation, from 2008 August 4 to 2012 May 
4, equating to a mission elapse time (MET) interval of 
239557417 to 357818882. In accordance with the PASS7 cri- 
teria, a zenith cut of 100° along with a rock-angle cut of 52° 
was applied to the data to remove any cosmic ray induced 
7-rays from the limb of the Earth's atmosphere. All 'source' 
class eventfl in a 15° radius of interest (Rol) centered on 
PKS 1510-089 were considered in the 0.1 GeV < E 1 < 300 
GeV energy range. 

Throughout this analysis, Fermi-LAT Science Tools 
version v9r27p1 were used in conjunction with instrument 
response functions (IRFs) p7source_v6. To investigate the 
7-ray properties of PKS 1510-089, I utilized the unbinned 
maximum likelihood estimator of the Fermi-LAT Science 
Tools' GTLIKE routine. GTLIKE allows us to fit the data with 
a series of both point and diffuse sources of 7-rays. The 
model used to calculate the likelihood of 7-ray emission from 



2 'Source' class events equates to EVENT_CLASS=2 in the PASS7 
data set. 



PKS 1510-089 was a combination of the most recent galactic, 
gal_2yearp7v6_v0.fits, and extragalactic, iso_p7v6source.txt, 
diffuse models, and all point sources within a 15° Rol cen- 
tered on PKS 1510-089. Each point source was modeled with 
a power-law spectrum of the form dN/dE — Ax (E / E )~ T . 
Both the photon index, V, and normalisation, A, of point 
sources within 10° of PKS 1510-089 were free to vary dur- 
ing the model fitting, while, for sources greater than 10° 
from PKS 1510-089, both T and A were fixed to the values 
publis hed in the Second Fermi Source Catalog (| Nolan et al.l 
(2Qlg)). 

3 7-RAY CHARACTERISTICS 
3.1 Light curve 

To investigate the temporal behavior of the 7-ray flux, the 
3.75 year data set was binned into daily temporal bins, with 
the GTLIKE routine applied to each bin separately. Only time 
intervals where the corresponding test statistic^, TS, was 
greater than 10 were considered, which equates to a signif- 
icance of ss 3cr. The resultant light curve can be seen in 
Figure 1, with statistical errors only. 



3 The test statistic, TS, is defined as twice the difference between 
the log-likelihood of two different models, 2[logL — logLo], where 
L and Lo are define d as the likelihood wh en the source is included 
or not respectively llMattox et al.l il996l) ). 



4 Anthony M. Brown 



(/) too 

>> 
n 
Q 



CD 
S3 

E 

3 



1 lH r i 




10" 5 

Flux (ph cm" 2 s 1 ) 



Figure 2. Upper panel: Duty Cycle of PKS 1510-089's 0.1 GeV 
< E-y < 300 GeV 7-ray flux during the first 3.75 years of Fermi- 
LAT observations, with the errors simply being the statistical 
error of the counts in each bin. Lower panel: Cumulative distri- 
bution function of the duty cycle, showing the number of days 
where the flux was below a given value F, normalised to the total 
number of days where the TS value was TS > 10. The cumu- 
lative distribution function clearly shows that the time spent by 
PKS 1510-089 at extreme flux levels of > 10" 5 photons cm" 2 
s" 1 constitutes a small fraction of the total observing period. 



Throughout the first 3.75 years of operation, PKS 1510- 
089 was one of the brightest AGN detected by the Fermi- 
LAT. Due to these high flux levels, requiring TS > 10 re- 
moves few daily bins from Figure 1. The variable nature 
of PKS 1510-089's 7-ray flux is clear to see, with the 3.75 
year light curve exhibiting short periods of flaring activity 
separated by extended periods of low fluxes on the order of 
(0.5-1) x 10" 6 photons cm" 2 s" 1 . PKS 1510-089's average 
daily 7-ray flux during the observing period is (1.39 ± 0.05) 
xlO -6 photons cm -2 s" 1 , though on 3 separate days, the 
flux exceeded 10" 5 photons cm -2 s" 1 at the peak of the flare 
period. During these flares, the 7-ray flux was at historical 
maxima for PKS 1510-089. 

To quantify the extent of the variability seen in Fig- 
ure [TJ the duty cycle of the 0.1 GeV < E 1 < 300 GeV 
7-ray flux during the 3.75 year observing period was de- 
rived. Defined as the fraction of time a source is at any 
given flux level, the duty cycle allows us to investigate if the 
obse rved flaring activ i ty is t y pical or aty p ical fo r the source 
(e.g. IVercellone et all 1(2004 ), lLott et all (|2012T l). Shown in 
Figure O the duty cycle utilised all TS > 10 flux measure- 
ments from the daily-binned light curve of Figure [T] The 
duty cycle's bins were defined by the minimum and max- 
imum daily flux values observed, with a total of 100 bins 
within the flux range. 

The 3.75 year duty cycle of PKS 1510-089 appears to 
have a well defined peak, where PKS 1510-089 spends the 
majority of time with a gradual decrease both above and 
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Figure 3. Upper panel: 6 hour binned light curve of the 0.1 
GeV < E-, < 300 GeV 7-ray flux from PKS 1510-089 during 
the 'flare 1' period, 55850 < MJD < 55858. Lower panel: 3 hour 
binned light curve of the 0.1 GeV < < 300 GeV 7-ray flux 
from PKS 1510-089 during the 'flare 1' period, 55850 < MJD 
< 55858. During this period, the flux exceeds 10 -5 photons cm -2 
s" 1 in the 3 hour binned light curve on five separate occasions. 
All temporal bins have a TS > 10. 



below this peak value, and an excess of events at the highest 
flux levels associated wit h the 3 days were the fl ux was above 
10~ 5 photons cm -2 s- MTavecchio et al.l l|201Ch constructed 
a Fermi-LAT duty cycle for PKS 1510-089 from 1.5 years of 
LAT observations, finding that the source was at a high flux 
level for approximately 1% of the total observing period, 
where 'high flux level' is defined as a flux that is a factor 
of 10 larger than the average flux level over the observing 
period. However, Tavecchio et al. noted that this percentage 
would most likely decrease over a larger observing period 
since they studied PKS 1510-089 specifically because it was 
in a high activity state during the observing period. This 
does indeed seem to be the case with the exceptional flare 
events observed during 55850 < MJD < 55880, which are 
^10 times the 3.75 year averaged flux, covering only ~ 0.2% 
of the total observing period reported here. 



3.2 Flux variability 

To search for rapid flux variability on sub-day timescales, 
the daily light curve shown in Figure [1] was utilised to select 
periods of high flux level. In particular, I concentrated on 
two flare periods which encompass the three days where the 
7-ray flux from PKS 1510-089 was at historical maximum. 
The first period, 'flare 1', spanning the period 55850 < MJD 
< 55858, saw the 0.1 GeV < E 1 < 300 GeV 7-ray flux 
from PKS 1510-089 peak at (1.53 ± 0.02) x 10" 5 photons 
cm" 2 s" 1 on MJD = 55854. The second period, 'flare 2', 
spanned the period 55866 < MJD < 55876, during which 
time, PKS 1510-089's daily flux twice exceeded 10 -5 photons 
cm" 2 s" 1 ; (1.28 ± 0.03) x 10" 5 photons cm" 2 s" 1 on MJD 
= 55868 and (1.23±0.08) x 10" 5 photons cm" 2 s" 1 on MJD 
= 55873. These high flux levels allowed us to re-analyse the 
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Table 1. Summary of quickest variability timescales events of PKS 1510-089 during the 3.75 year period, which 
are less than 3 hours and have a significance of at least 5<r. The times, T start and Tstop, are in MJD, with the 
fluxes in units of 10 — 6 photons cm -2 s _1 . The observed characteristic timescale r, from Equation 1, is converted 
to the intrinsic variability timescale, Tint, with Tint = t(1 + z)~ 1 . The uncertainty in the variability timescale was 
calculated by propagating the uncertainty in the flux and time values through Equation 1. The last column indicates 
whether the variability event is an increase (rise) or decrease (decay) in the flux. 
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observations with GTLIKE, in 6 and 3 hour bins and still 
satisfy the TS > 10 criteria for the majority of temporal 
bins. The resultant light curves can be seen in Figure for 
flare 1 and Figure [4] for flare 2. 

As can be seen in both Figure [3] and Figure [4j the 3 
hour binned light curve reveals a large amount of variability 
that was masked, or washed out, by the daily bins utilised in 
Figure [l] Furthermore, the smaller temporal bins of Figure 
[3] and Figure [4] reveal that the flux exceeded 10 -5 photons 
cm -2 s _1 during 17 separate 3 hour intervals, compared 
to 3 daily periods revealed in Figure [1] with the 0.1 GeV 
< E 1 < 300 GeV 7-ray flux from PKS 1510-089 peaking at 
(1.99±0.04)xl0~ 5 photons cm" 2 s"\ equating to ~ 3 x 10 48 
ergs s _1 , on MJD= 55872.937 ± 0.063. 

To characterise the timescales of the observed flux vari- 
ability, the time taken for the flux to increase or decrease 
by a factor of 2 was evaluated. Referred to as the dou- 
bling/halving timescale depending upon whether the flux 
is increasing or decreasing, this timescale is defined by: 

f(i) = F((„)x2 (T_1(t " w) (1) 

where r is the characteristic doubling or halving 
timescale and F(t) and F(t ) are the fluxes at time t and 
t respectively. The above equation was applied to consecu- 
tive flux measurements that satisfy the TS > 10 criteria, as 
shown in Figure [T] Figure [3] and Figure [4l in a systematic 



search for the quickest variability timescale of PKS 1510- 
089's 0.1 GeV < E 1 < 300 GeV 7-ray flux. No selection 
criteria is applied to the flux level and its associated error; 
rather, only resultant variability timescales that have a value 
which is greater than 5 standard deviations from zero were 
considered. A summary of intrinsic variability timescales 
with I Ti n t \< 3 hours, which have a significance of at least 
5a, are given in Table 1. 

This approach has two important caveats: firstly, for 
consecutive flux values where the difference is less than a 
factor of two, the variability timescale calculated is essen- 
tially a prediction of how long it would take for a double or 
halving of flux to occur if the observed change in flux contin- 
ued at its current rate. Secondly, the variability timescales 
calculated assumes that the flux increase or decrease is con- 
stant throughout the 3 hour bin, and as such, should be con- 
sidered as an upper limit. Probing timescales smaller than 
those observed is primarily limited by the 3 hour period that 
Fermi-LAT takes to scan the entire sky and the amount of 
time a source is in the bore sight of the LAT instrument 
where the sensitivity is the greatest. 

As one would expect, the 3 hour binned light curves 
of Figure [3] & [4] reveal the most rapid variability. In par- 
ticular, the shortest doubling rise time, T r i 3e , was found to 
be 1.30 ± 0.12 hours at MJD= 55875.437 ± 0.063, and the 
shortest halving decay time, Td ecay , was 1.21 ± 0.15 hours 
at MJD= 55852.187 ± 0.063. It is worth highlighting that 
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R, with a Doppler factoiQ S, is related to the 7-ray variability 
timescale, t var , by: 



i66 55868 55870 55872 55874 55876 

MJD 

Figure 4. Upper panel: 6 hour binned light curve of the 0.1 
GeV < E-, < 300 GeV 7-ray flux from PKS 1510-089 during 
the 'flare 2' period, 55866 < MJD < 55876. Lower panel: 3 hour 
binned light curve of the 0.1 GeV < E 1 < 300 GeV 7-ray flux 
from PKS 1510-089 during the 'flare 2' period, 55866 < MJD 
< 55876. During this period, the flux reached a historic maximum 
for PKS 1510-089 of (1.99 ±0.04) X 10" 5 photons cm" 2 s" 1 on 
MJD= 55872.937. All temporal bins have a TS > 10. 



there are several instances of variability on even smaller 
timescales, for example on MJD= 55852.437 ± 0.063 the in- 
trinsic variability timescale was —0.91 ±0.66 hours, however 
with large uncertainties in the variability timescale, these 

events do not satisfy the > 5a criterion. 

A similar study was undertaken by iFoschini et al.l 

|201ll ) for the FSRQs 3C 454.3, 3C 273 and PKS B1222±216 
utilising ~ 2 years of LAT observations. This study revealed 
intrinsic flux variability in the 2 — 3 hour range, but only 
at the 3<r level of certainty. Furthermore, Foschini et al. 
adopted a conservative approach, and only calculated the 
variability timescale between subsequent flux measurements 
with a 3a difference in magnitude. If the same approach 
is applied to the variability events in Table 1, the shortest 
rise and decay timescales are still found to be 1.30 ± 0.12 
and 1.21 ± 0.15 hours respectively. Indeed, even requiring 
a 5cr difference in consecutive flux values to calculate the 
associated doubling/halving timescale, result in the same 
timescale for the quickest decay in flux, and only a slightly 
longer rise time of 1.36 ±0.13 hours. As such, the 1.21 hour 
decay-time and 1.3 hour rise-time timescales observed in 
this study are the quickest flux variability observed from 
the FSRQ subclass of AGN in the MeV-GeV energy range. 
The implications for the size of this variability is discussed 
in §4. 



R < ct var 5(l + zy 



(2) 



where z is the redshift of the source. For the intrinsic 
variability timescales outlined in Table 1, the size of the 7- 
ray emission region for the observed flares can be constrained 
to be R5 _1 < (0.9—2.3) x 10 14 cm. For comparison, adoptin g 
a SMBH mass of 5.4 x 10 8 M from lAbdo et alj (|2010bl) . 
the Schwarzschild radius of PKS 1510-089 is ~ 1.6 x 10 14 
cm. 

Radio observations of PKS 1510-089 during February 
to June 2010 have shown that structures w ithin its rela- 
tivisti c jet have Doppler factors of S = 47 (|Kadota et al.l 
(2012)), while observations d uring 1998 to 2 001 s how sim- 
iliar values of 5 = 34 - 42 jjorstad et alj (|2005l ) ). If one 
assumes these rather extreme values for the Doppler factor 
are a good representation of the average velocity structure 
within PKS 1510-089's relativistic jet, the emission regions 
of the flares reported here are found to be less than 3 x 10~ 3 
parsecs in size. This suggests that only a very small portion 
of the relativistic jet is responsible for the observed 7-ray 
flares. 

While these rapid variability timescales allow us to gain 
insights into the physical properties of the 7-ray emission 
region, the extreme nature of the variability does little to 
constrain the location of the emission region. Small structure 
within the relativistic jet is not unique to just the inner 
regions of the AGN; small over-densities of plasma can occur 
in the jet at both the sub-par s ec, BLR, and the parsec-scale , 
MT region dNishikawa et alj (l2003t). iPerucho et al.1 (|2006h . 
iBromberg fc Levinsonl (|2009l ), iKohler et alj (|2012r i). 

External to the jet however, the environment of the 
BLR and MT are significantly different, with the energy dis- 
tribution of the BLR photon field peaking at UV energies, 
while in the MT it peaks at IR energies. This difference in 
photon energies, and the subsequent cooling effect they have 
on the energy distribution of the particle population within 
the emission region, can be utilised to constrain the origin 
of t he 7-rays. 

iTavecchio et al.l (|2010l ) exploited this difference to con- 
strain the cooling timescales for a relativistic electron pop- 
ulation located within the BLR and within the MT region. 
Utilising Equation 3 below, Tavecchio et al. calculated the 
cooling timescale for an electron producing a 100 MeV 7- 
ray photon, through the Inverse-Compton (IC) process, to 



be t° 



800 seconds in the BLR and t° 



12, 000 sec- 



onds in the MT. Furthermore, for a 5 keV X-ray photon 
produced via the IC process, Tavecchio et al. calculated an 



electron cooling timescale of t 



and t° 



20 days in the M 



31 hours in the BLR 



4 ORIGIN OF THE 7-RAY EMISSION 
4.1 Variability timescale 

Taking the Doppler factor of the relativistic jet into consid- 
eration, causality implies that the size of an emission region, 



4 8 = (r(l — /3cos)) — 1 where V is the bulk Lorentz factor of the 
jet, /3 = v/c and 6 is the angle to the line of sight. 

5 Given that X-ray cooling timescales from blazar s have been ob- 
served on the sub-hou r timescale, (fo r example, see I Foschini et ahl 
i2006h for FSRQs and lBrownl J2006T) for BL Lac objects), the ap- 
plication of this approach to observed X-ray flux variability is 
questionable. 
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55850 55851 55852 55853 55854 55855 55856 55857 55858 

HMD 

Figure 6. Upper panel: the 0.1 GeV < E~ < 1 GeV 7-ray flux 
from PKS 1510-089 during the 'flare 1' period, 55850 < MJD 
< 55858, binned in 6 hour periods. Lower panel: the 1 GeV < 
E-y <300 GeV 7-ray flux from PKS 1510-089 during the 'flare 
1' period, 55850 < MJD < 55858, binned in 6 hour periods. All 
temporal bins have a TS > 10. 
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Figure 7. Upper panel: the 0.1 GeV < E-y < 1 GeV 7-ray flux 
from PKS 1510-089 during the 'flare 2' period, 55866 < MJD 
< 55876, binned in 6 hour periods. Lower panel: the 1 GeV < 
E 1 <300 GeV 7-ray flux from PKS 1510-089 during the 'flare 
2' period, 55866 < MJD < 55876 binned in 6 hour periods. All 
temporal bins have a TS > 10. 



ob3 _ 3m e c(l + 2) 

coot ~ 4a T U'yS (> 

The 3 — 6 hour cooling timescales that Tavecchio et 
al. observed did not allow them to distinguish between a 
BLR or a MT origin for the e mission region. H o wever util- 
ising the approach outlined in iTavecchio et alj (|20T0t ), the 
Tdecay = 1.21 ± 0.15 hour cooling timescale discovered in 
this study would suggest that the emission region, for the 
quickest flaring event, is located within the BLRq However, 
it is important to highlight that this approach assumes that 
the IC interaction is occuring in the Thomson regime, which 
as discussed in §4.2, is mostly likely not the case. As such, 
this approach might not give a realistic cooling timescale, 
which would bring the conclusion about the location of the 
emission region into question. 

4.2 Energy dependent cooling timescales 

Tavecchio et al. assumed that the IC scattering was occur- 
ing in the Thomson regime, with *ye « 1 being true for all 
electron-photon interactions, where 7 is the Lorentz factor 
for the electron and e is the energy of the photon. Dotson 
et al. (2012a,b) on the other hand, found that the IC scat- 
tering occurred in the Klein-Nishina regime when the emis- 
sion region was located in the BLR, while it occurred in 
the Thomson regime for emission regions located within the 
MT region of the AGN. This difference manifests itself as 
an energy dependent cooling timescale for emission regions 

6 It should be noted that this conclusion does not apply to all 
variability timescales reported in Table 1 since, for events with 
Tdecay > 2.5 hours, it is not possible to differentiate between 
the BLR and MT cooling timescales of 800 and 12,000 seconds 
respectively. 



embedded in the MT region and a (quasi-) energy indepen- 
dent cooling timescale for emission regions embedded in the 
BLR. Specifica lly, utilising the BLR and MT photon energy 
densities from iGhisellini fe Tavecchid ([2009), Dobson et al. 
calculated a < 1 hour difference in the flux halving timescale 
for 200 MeV photons and 20 GeV photons for an emission 
region located in the BLR, but a ~ 10 hour difference in 
flux halving timescales for the two photon en ergies from an 
emissio n region embedded in the MT region |Dotson et al.1 
(|2012bh ). 

More generally speaking, this difference would result 
in a time lag between the cooling of the MeV and GeV 
components of a 7-ray flare. Conversely, if present, a time 
lag can be u sed to constra i n the e nergy density of the MT 
photon field (|Dotson et al.1 l|2012ah l. 

To investigate the possibility of energy dependent cool- 
ing timescales, the daily light curve of Figure 1, along with 
the 6 hour binned light curves of Figure 3 and Figure 4, were 
re-analysed as per the procedure outlined in §2, but for pho- 
tons in two distinct energy groups; low-energy, 0.1 — 1 GeV, 
and high-energy, 1 — 300 GeV. Due to the limited number 
of events above 1 GeV and the requirement that each tem- 
poral bin has TS > 10, this re-analysis was not applied to 
the 3 hour binned light curve. The resultant daily binned 
light curves for both low and high energy fluxes, can be seen 
in Figure [5] The low and high energy light curve, for both 
flare 1 and flare 2, binned in 6 hour intervals, are shown in 

Figure [5] and Figure [7] respectively. 

Discrete correlation functions, fDCF; lEdelson fc Krolikl 

(1988)), were applied to the low and high energy light curves 
of Figure HJ] and Figure [7] to search for the presence of a time 
lag in the flux level changes between the energy bands. An 
important characteristic of the DCF procedure is that it 
allows us to accommodate for differences in the sampling 
rates of the two light curves associated with the TS > 10 
criterion. The DCFs represent a more robust approach to 
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Figure 5. C/pper panel: Lightcurve of the 0.1 GeV < £ 7 < 1 GeV 7-ray flux from PKS 1510-089, during the period 2008 August 4 (MJD 54682) 
to 2012 May 4 (MJD 56051), binned in daily periods. Lower panel: Lightcurve of the 1 GeV < E 1 <300 GeV 7-ray flux from PKS 1510-089, 
during the period 2008 August 4 (MJD 54682) to 2012 May 4 (MJD 56051), binned in daily periods. All temporal bins have a TS > 10. 



searching for energy (in) dependent cooling timescales than 
simply applying Equation 1 to the individual light curves 
since the latter is limited by the statistics and associated er- 
ror of the high energy light curve. Indeed, applying Equation 
1 to the daily binned light curves of Figure [5] does not find 
any 1 — 300 GeV flux cooling timescales at a 5a level of confi- 
dence. Furthermore, if the application of Equation 1 reveals 
the presence of significant variability timescales in both the 
0.1 — 1 and 1 — 300 GeV light curves, a-priori knowledge of 
the magnitude of the time lag is needed to determine if the 
variability timescales from the individual energy ranges are 
related. 

The resultant DCFs for the flare 1 and flare 2 events 
are shown in Figure [8] and Figure [^respectively, binned in 6 
hour intervals. A positive time lag, i(o.i-iGeV) — i(i-300GeV)> 
between the two light curves implies that the 0.1 — 1 GeV 
flux is delayed with respect to the 1 — 300 GeV flux. However, 
it is important to note that this delay applies to both the 
flux increases and flux decreases. To interpret any observed 
lag in the DCF as evidence for energy dependent cooling 
timescales, one has to assume that there the flux increase in 
both energy bands occurs at the same time. 

The DCF for the flare 2 event exhibits a clear peak. A 
gaussian fit to this peak indicates that this temporal lag is 
—6.13 ± 2.83 hourfl However the location of the peak sug- 
gests that changes in the 0.1 — 1 GeV flux in fact proceed 



any changes in the 1 — 300 GeV flux. A closer inspection of 
the light curves in Figure [7] finds that the 0.1 — 1 GeV flux 
peaks approximately 6—12 hours before the 1 — 300 GeV 
flux. As such, it would appear that, if an energy dependent 
cooling timescale is present, it is masked by the energy de- 
pendent rise time of the flare event. 

The DCF for flare 1 on the other hand, appears to have 
a peak at a time-lag of 0. This would suggest that there is 
no energy dependence to the cooling timescale, and therefore 
point towards a BLR-origin for the 7-ray emission. However 
there is a large amount of uncertainty in the DCF for neg- 
ative time-lags, with additional peaks at ~ —2 and ~ —2.5, 
and thus it is not possible to draw strong conclusions about 
the DCF peak at and the DCF for flare 1 in general. Fur- 
thermore, the DCF for the flare 2 event suggests that the 
DCF could be dominated by effects not related to the elec- 
tron cooling timescales. As such, this implies that for these 
flare events, little is gained from the DCF approach outlined 
here with regards to determine the location of the 7-ray 
emission region. 



4.3 Photon-photon pair-production 

Another important difference between the two possible loca- 
tions for the 7-ray emission regions is the role that photon- 
photon pair-production, (77 — ► e + e~) plays in attenuat- 



7 It is worth noting that this peak feature is still present if the 
DCF is rebinned into 12 hour intervals, while the 12 hour binned 



DCF for flare 1 still has no well defined peak, and large errors for 
negative time lags. 
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Figure 8. Discrete correlation function calculated between the 
0.1-1 GeV and the 1 - 300 GeV light curves of Figure \E\ for the 
flare 1 period, binned in 6 hour intervals. The shading indicates 
the error of the DCF. The DCF indicates that no time lag is 
present between the two light curves during this flare period. 
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Figure 10. The difference in TS values, defined as TS^ — 
TS pome) — l aw , where x is the TS value for either the broken power- 
law, a broken power-law with an exponential cut-off and a log- 
parabola fits, throughout the 'flare 1' period. To improve statis- 
tics, especially at the high-energy tail of the spectrum, the differ- 
ence in model TS values was calculated on daily intervals of obser- 
vations. There is a TS > 25 difference between a power-law and 
a log-parabola description of PKS 1510-089 during MJD= 55853. 
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Figure 9. Discrete correlation function, calculated between the 
0.1-1 GeV and the 1 - 300 GeV light curves of Figure [7] for the 
flare 2 period, binned in 6 hour intervals. The shading indicates 
the error of the DCF. A positive time lag implies that the 0.1 — 1 
GeV flux is delayed with respect to the 1 — 300 GeV flux. A 
gaussian fit to the DCF peak finds a time lag of —6.13 hours to 
be present. 



ing the emitted 7-ray flux and the subsequently imprint 
that the process leaves o n the observed 7 -ray spectrum 
l|Donea fc Protherod <|2003h . lLui &~Bail (2006)). Simply put, 
as a result of photon-photon pair-production, the BLR of 
FSRQs is opaque to 7-rays above ~ 20 GeV in energy, while 
the MT region is not. As a result, if the 7-ray region is lo- 
cated in the BLR, one would expect to see a cut-off in the 
7-ray spectrum attributed to pair production, while a 7-ray 
spectrum originating from the MT would not have such a 
feature. 

To search for a spectral cut-off feature, the flare 1 and 



flare 2 light curves were again re-analysed in daily bins, with 
PKS 1510-089 modeled by 3 additional models to the power- 
law used originally in §2. To improve statistics, especially 
at the high-energy tail of the spectrum, the observations 
were re-analysed in daily intervals. The additional models 
utilised were a broken power-law, a broken power-law with 
an exponential cut-off and a log-parabola. All other point 
and diffuse sources in the Rol were modeled as before (see 
§2). The TS value for each model, in each daily bin, was 
compared to that of a simply power-law, to investigate if 
there is a significant deviation from a power-law spectrum 
during the flaring event. The difference in TS is defined as 
TS;e— TSpowei — law, where x is the TS value for either the 
broken power-law, broken power-law with an exponential 
cut-off or log-parabola fits. A large positive difference in TS 
values indicates that the 7-ray spectrum is better described 
by a function with a cut-off feature then a simple power-law 
distribution, and as such, suggests the presence of a spectral 
cut-off. An important caveat of this approach is that the 
presence of a spectral cut-off does not automatically imply 
a BLR origin for the 7-ray flux; for example a cut-off in 
the 7-ray spectrum can also occur if there is a cut-off in the 
energy distribution of the particle population responsible for 
the 7-ray emission. 

The daily differences in TS values during the flare 1 and 
flare 2 events can be seen in Figurc[l0]and Figure [TT] respec- 
tively. Flare 1 exhibits a large TS difference, ATS > 25, 
between a power-law and a log-parabola description of PKS 
1510-089's 0.1 GeV < E < 300 GeV 7-ray spectrum during 
MJD= 55853, a day before flare l's maximum flux. How- 
ever, on MJD= 55850, at the start of flare 1, there is little 
difference in the TS value from either a power-law or log- 
parabola description of PKS 1510-089's spectrum. Likewise, 
towards the end of the flare 1 event, there is also little differ- 
ence in the TS values of the different spectral shapes. This 
variation in the TS value for the log-parabola description 
of the spectrum, relative to the power-law fit, suggests that 
throughout the flare 1 event, there is a change in the shape 
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Figure 11. The difference in TS values, defined as TSa; — 
TS power -i aw , where x is the TS value for either the broken 
power-law, a broken power-law with an exponential cut-off and 
a log-parabola fits, throughout the 'flare 2' period. To improve 
statistics, especially at the high-energy tail of the spectrum, the 
difference in model TS values was calculated on daily intervals of 
observations. No significant difference in TS value is seen through- 
out the second flare. 



of the 7-ray spectrum. What is more, the preference of the 
log-parabola fit over the power-law fit implies the presence 
of a cut-off in the spectrum. 

To investigate this spectral cut-off, the 7-ray spectrum 
during MJD= 55850 and MJD= 55853 were fitted with both 
a power-law and log-parabola function, and the reduced x 2 
of the function fit was calculated. The spectra were obtained 
by applying GTLIKE separately to f logarithmic energy bins 
in the 100 MeV to 102.4 GeV energy range. Only energy bins 
with a GTLIKE TS value greater than 10 were considered in 
the function fit. The resultant spectrum, power-law and log- 
parabola fits for MJD= 55850 and MJD= 55853 can be seen 
in Figure [121 and Figure [131 respectively. 

From Figure 1121 it can be seen that at the start of the 
flare 1 event, when there is no difference in the TS values, 
the spectrum is well described by a power-law, with no in- 
dication of a spectral cut-off presen10. However the 7-ray 
spectrum during MJD= 55853, seen in Figure 1 1 31 exhibits 
a clear deviation from a simple power-law distribution at 
the high-energy tail of the spectrum. The deficiency of the 
power-law is confirmed by the reduced \ 2 °f the power-law 
fit, which is 1.885, while the reduced \ 2 °f the log-parabola 
fit is 0.282. The reduced \ 2 °f the log-parabola fit thus con- 
firms the presence of a cut-off in the 7-ray spectrum for the 
flare 1 event, and therefore suggests that the flare 1 event 
originates from within the BLR. 

For flare 2, on the other hand, Figure [TT1 indicates that 
there is no significant deviation from a power-law description 
throughout the flare event, with the greatest difference in TS 
values being ~ 10. What is more, for a large number of days, 
the TS value for the power-law fit is greater than that of the 
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Figure 12. The spectrum of PKS 1510-089 during MJD= 55850. 
The solid green line indicates the best power-law fit to the spec- 
trum, while the dash red line indicates the best log-parabola fit to 
the spectrum. The reduced x 2 of the power-law fit is 0.112 while 
the reduced x 2 of the log-parabola fit is 0.155. The GTLIKE fit to 
the 3.2 — 6.4 GeV energy band has a TS value less than 10, and 
as such, was replaced by a 95 % confidence upper limit, however 
this limit was not included in the reduced x 2 calculation for the 
function fit. 
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Figure 13. The spectrum of PKS 1510-089 during MJD= 55853. 
The solid green line indicates the best power-law fit to the spec- 
trum, while the dash red line indicates the best log-parabola fit to 
the spectrum. The reduced x 2 of the power-law fit is 1.885 while 
the reduced x 2 of the log-parabola fit is 0.282. 



spectral models with cut-offs. As such, the lack of a spectral 
cut-off suggests an MT origin for the 7-ray emission. 



8 It is important to note that the log-parabola function, which is 
defined by |g = A X £("-/3xio S i0(E)) _ also describes the spec _ 

trum well, simply because the co-efficient /3 is close to zero, thus 
rendering the log-parabola function, a power-law function. 



5 MULTIPLE EMISSION REGIONS 

Taken in isolation, the presence or lack of a spectral cut-off 
or temporal lag in the discrete correlation functions, do not 
provide any strong evidence as to the location of the emis- 
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sion region. However, the combination of the two properties 
allows for a stronger argument as to the possible location. 

Flare 1: At its peak flux, the 0.1 — 300 GeV spectrum 
of flare 1 (55850 < MJD < 55858) showed significant devi- 
ation, greater than TS > 25, from a power-law description, 
indicating the presence of a spectral cut-off. The presence 
of this cut-off was confirmed by the reduced \ 2 values for 
the power-law and log-parabola fits to the 0.1 — 300 GeV 
spectrum on MJD= 55853, with the power-law fit having 
a reduced x 2 = l-885 and the log-parabola fits having a re- 
duced x 2 =0.282. Due to photon-photon pair-production, the 
spectral cut-off suggests that the 7-ray emission region as- 
sociated with this flare is located within the BLR, on the 
sub-parsec scale from the central SMBH. While a peak at a 
time-lag of in the 0.1-1 GeV versus 1 - 300 GeV DCF 
for flare 1 would suggest that there is no energy depen- 
dence to the cooling timescale, and thus a BLR-origin for 
the 7-ray emission, the uncertainty in the DCF for negative 
time-lags is too large to draw strong conclusions. Further- 
more, the DCF for the flare 2 event suggests that the DCF 
is dominated by effects not related to the electron cooling 
timescales and as such, implies that little is gained from the 
DCF approach outlined here with regards to determining 
the location of the 7-ray emission region. 

Flare 2: During the the flare 2 period (55866 < MJD 
< 55876) the daily flux twice exceeded 10 -5 photons cm -2 
s _1 . However there was no significant deviation from a 
power- law distribution at any point during the flare event. 
As such, without invoking the presence of axion-like parti- 
clefl, the lack of spectral cut-off points towards a MT loca- 
tion for the 7-ray emission region associated with this flare. 
The 0.1-1 GeV versus 1 - 300 GeV DCF for flare 2 does not 
reveal any evidence for an energy dependence to the cooling 
timescales, though this may to be due to the fact that the 
DCF is dominated by an apparent energy dependence in the 
flare rise time. 

While there is ~ 8 days between the two flaring 
episodes, there is evidence indicating that the two flare 
events are spatially separated by > 1 parsec, with flare 1 be- 
ing of BLR origin and flare 2 being of MT origin. Thus a nat- 
ural conclusion to draw from these studies is that there are 
multiple, simultaneously active, 7-ray emission zones along 
the relativistic jet in both the BLR and MT, capable of 
emitting 0.1 — 300 GeV 7-rays. 

Multi-zone emission models have often been used to ex- 
plain the broad-b and spec t ral en e rgy distribution o f blaza rs 
(for example, see iBrownl | |2006l ). iNalewaiko et all |2012l )'). 
or to explain the 7 -ray emission from misaligned AGN 
|Lenain et alj |200ct )L However, for these multi-zone mod- 
els, the 7-ray emission origin is either in the BLR or the 
MT, never both. 

From a multi-wavelength (MWL) campaign of PKS 
1510-089, during the first 11 months of Fermi-LAT oper- 
ation, Abdo et al. conclude d that the 7-ray em ission origi- 
nated from within the BLR |Abdo et al ] (|2010bl )1; this con- 

9 If the flare 2 event was in fact located within the BLR, 
axion-like particles could possibly explain the absence of a spec - 
tral cut-off (e.g see iTavecchio et al.1 J2012h . iHorns et al.1 1 120121) . 
iRoncadelli et aI.N2012rO . Without axion-like particles, the combi- 
nation of BLR origin and lack of a spectral cut-off would require a 
sudden drop in the local energy density of the BLR photon field. 
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Figure 14. The black crosses represent hardness ratio, defined 
as the 1 - 300 GeV flux divided by the 0.1 - IGeV flux, of PKS 
1510-089 over the entire 3.75 year period. The red filled triangles 
represent the arrival times of all photons > 20 GeV in energy 
within 0.1 degrees of PKS 150-089. The scale for the hardness 
ratio is on the left, while the photon energy scale is on the right. 
There does not appear to be any preference for the highest energy 
photons to arrive during hard spectral states. For comparison, 
the MJD when VHE emission was detected by the ground-based 
HESS and MAGIC telescopes are indicated by the vertical arrows. 

elusion primarily being driven by the presence of a cut-off 
in the 7-ray spectrum. However, both Marscher et al. and 
Orienti et al. have concluded from their respective MWL 
campaigns of PKS 15 10-089 that the 7-ray o riginates from 
outsid e BLR region (jMarscher et al.l (|2010| ). lOrienti et al.l 
I|2012r 0. with their conclusions primarily based on the flaring 
events being observed simultaneously at 7-ray, optical and 
radio wavelengths. The presence of multiple 7-ray emission 
regions along the jet allows us to reconcile these seemingly 
contradictory conclusions, with PKS 1510-089 being able to 
produce 7-ray flares from regions both within the BLR and 
the MT regions of it's relativistic jet. 



6 VHE EMISSION 

Finally I briefly turn my attention to emission of Very High 
Energy (VHE) 7-rays from PKS 1510-089. The detection 
of VHE 7-ray photons from FSRQs is difficult to accom- 
modate in a pure BLR-origin model for the location of the 
7-ray emission region. VHE 7-ray emission from PKS 1510- 
089 was originally detected by the HESS telescope array 
in 2010 and ha s been r e cently confirmed by the MAGIC 
telescope array (|Wagnerl (|20ld h [Cortina! (|2012h '). A gtlike 
likelihood analysis of all > 50 GeV photons from PKS 1510- 
089 detected by .Fermi-LAT, utilizing the same point source 
and diffuse model as that of §2, finds PKS 1510-089 to be 
a source of VHE 7-rays at the ~ 4.6<r confidence level. Fur- 
thermore, a gtlike likelihood analysis finds PKS 1510-089 
to be a source of > 20 GeV 7-rays at the ~ 13.5<r confidence 
level. 

A closer inspection of the individual photon events re- 
veals that the highest energy photon detected from PKS 
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1510-089 has an energy of 70 GeV, with a total of 11 
events above 2 GeV withi n 0.1 degrees of PKS 1510 -08 jB- 
iKataoka et"aH l|2010h and H rown fc Adamj (|201lh found 
that it is the >GeV 7-ray flux and spectral shape that is 
important when triggering ground-based VHE 7-ray obser- 
vations, with a higher >GeV flux or, harder 7-ray spectrum, 
more likely to be associated with the emission of VHE 7-ray 
photons. To investigate if this applies to PKS 1510-089, the 
arrival times of the > 20 GeV photons was compared to 
the hardness ratio of the light curves of Figure [5] The hard- 
ness ratio was defined as the 1 — 300 GeV flux divided by 
the 0.1 — 1 GeV flux. The light curve of the hardness ratio, 
binned in daily intervals, along with the arrival times of the 
> 20 GeV photons, can be seen in Figure [14] 

There is no evidence in Figure Q3] of a correlation be- 
tween the arrival of > 20 GeV photons and a spectral hard- 
ening as shown by a higher hardness ratio. Furthermore, 
there does not appear to be any obvious temporal cluster- 
ing of > 20 GeV photons and the hardness ratio only ex- 
ceeds 0.2 on one occasion. Interestingly, there are no > 20 
GeV photons around the time where there is a spike in the 
hardness ratio of ~ 0.6. 

The MAGIC detection of PKS 1510-089 resulted from 
10 hours of observations taken during the period 2012 Febru- 
ary 3 to 2012 February 20 (55960 < MJD < 55980), either 
side of a full moon period l|Cortinal (|2012l )). From Figure [14] 
we see that there is no obvious increase in the hardness ratio 
during this period, and no > 20 GeV photons were detected 
either. Likewise, during the HESS observations, March 2010, 
only one > 20 GeV photon was detected, and the hardness 
ratio exhibits a slight decreasing trend. As such, it would 
appear that for PKS 1510-089, there is no obvious 0.1 — 300 
GeV flux trend that leads to the emission of VHE 7-rays. 
This can be interpreted as further evidence of multiple emis- 
sion regions along the jet being simultaneously active, since 
such a scenario would lead to complex variability patterns, 
resulting in an apparent random variation of the hardness 
ratio, that would mask an underlying pattern of VHE emis- 
sion occuring when the >GeV 7-ray flux increases. 



7 CONCLUSIONS 

An in-depth study of the 0.1 — 300 GeV 7-ray properties of 
PKS 1510-089 utilising Fermi-LAT observations from Au- 
gust 2008 to May 2012 has been reported. During this pe- 
riod, PKS 1510-089 has persistently been one of the bright- 
est and most variable AGN observed by the LAT detector. 
Several large 7-ray flares were observed, where the daily 
0.1 — 300 GeV 7-ray flux exceeded 10 -5 photons cm" 2 s — 1 . 
The photon statistics associated with these large flares al- 
lowed the daily-binned light curve to be re-analysed in both 
6 and 3 hour intervals and still satisfy a TS > 10 criteria for 
the vast majority of bins. The 3 hour binned light curve re- 
vealed the presence of doubling/halving timescales as small 
as T r i Se = 1.30 ± 0.12 hours and Tdecay = 1-21 ± 0.15 hours. 

The 10~ 5 photons cm' 2 s" 1 flare periods were studied 
in more detail in an attempt to find evidence for the loca- 
tion of PKS 1510-089's 7-ray emission region. In particular, 

10 0.1 degrees is larger than the 95% containment angle for pho- 
tons > 20 GeV. 



two approaches were used: (i) searching for an energy de- 
pendence to the cooling timescales, and (ii) searching for 
evidence of a spectral cut-off. 

The flare 1 event (55850 < MJD < 55858) possessed a 
spectral cut-off at its peak flux. However, the uncertainty in 
the DCF did not allow any clear conclusions on whether an 
energy dependency in the variability timescales is present. 
While there is an apparent peak in the DCF for flare 1 at a 
time-lag of 0, possibly hinting at a BLR origin for the flare, 
the large uncertainty in the DCF at negative time-lag values 
does not allow for any strong conclusions to be drawn from 
the DCF with regards to the origin of the 7-ray emission. 

The flare 2 event (55866 < MJD < 55876) showed no 
evidence for a spectral cut-off, possibly hinting at a MT 
origin of the 7-ray emission. If this is so, one would expect 
to see an energy dependency to the cooling timescale. While 
no evidence for this could be seen in the DCF of the 0.1 — 1 
GeV and 1 — 300 GeV flux levels, this may to be due to 
the fact that the DCF is dominated by an apparent energy 
dependence in the flare rise time. 

These seemingly contradictory origins for the observed 
7-ray flares can be reconciled if there are multiple, simulta- 
neously active 7-ray emission zones along PKS 1510-089's 
relativistic jet, capable of emitting 0.1 — 300 GeV 7-rays. 
Since the inner regions of FSRQs are opaque to 7-rays above 
10 — 20 GeV in energy, this multi-zone model can explain 
the existence of VHE 7-rays for PKS 1510-089, with emis- 
sion regions in the MT producing the observed VHE flux, 
while those inside the BLR contributing to the 0.1 — 10 GeV 
flux. While this study has unvealed evidence that both the 
BLR and MT location models are correct, further studies 
are needed to see if this property is unique to PKS 1510-089 
or a characteristic of FSRQs in general. 
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